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1. Introduction 

The Helton- Vinnikov Theorem [16 gives a geometric characterization of two- 
dimensional spectrahedra. They are precisely the subsets of M 2 that are bounded by 
rigidly convex algebraic curves, here called Helton-Vinnikov curves. These curves 
are cut out by hyperbolic polynomials in three variables, as discussed in [TS] . This 
theorem is a refinement of a result from classical algebraic geometry which states 
that every homogeneous polynomial in three variables can be written as 

f{x,y,z) = det(Ax + By + Cz) (1) 

where A 7 B and C are symmetric matrices. Here the coefficients of / and the matrix 
entries are complex numbers. When the coefficients of / are real then it is desirable 
to find A, B and C with real entries. The representations relevant for spectrahedra 
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are the real definite representations, which means that the linear span of the real 
matrices A, B and C contain a positive definite matrix. Such a representation is 
possible if and only if the corresponding curve {(x : y : z) € : f{x, y, z) = 0} is 
rigidly convex. This condition means that the curve has the maximal number of 
nested ovals, namely, there are d/2 resp. (d— l)/2 nested ovals when the degree d 
of / is even resp. odd. The innermost oval bounds a spectrahedron. 

Two linear matrix representations Ax + By + Cz and A'x + B'y + C'z of 
the same plane curve are said to be equivalent if they lie in the same orbit under 
conjugation, i.e. if there exists an invertible complex matrix U that satisfies 

U ■ {Ax + By + Cz) ■ U T = A'x + B'y + C'z. 

We call an equivalence class of complex representations real (resp. real definite) 
if it contains a real (resp. real definite) representative. Deciding whether a given 
complex representation is equivalent to a real or real definite one is rather difficult. 

We shall see that the number of equivalence classes of complex representations 
([!]) is finite, and, for smooth curves, the precise number is known (Thm. 2.1 1. Using 
more general results of Vinnikov [53] , we also derive the number of real and real 
definite equivalence classes. If a Helton- Vinnikov curve is smooth then the number 
of real definite equivalence classes of equals 2 9 , where g = C*^ 1 ) is the genus. 

This paper concerns the computational problem of constructing one represen- 
tative from each equivalence class for a given polynomial f(x, y, z). As a warm-up 
example, consider the following elliptic curve in Weierstrass normal form: 

f(x,y,z) = (x + ay)(x + by)(x + cy) - xz 2 . 

Here a, b, c are distinct non-zero reals. This cubic has precisely three inequivalent 
linear symmetric determinantal representations over C, given by the matrices 





x + cy 



x + ay 








x + by z 




x + cy 



All three matrices are non-real if o, b and c have the same sign, and otherwise 
two of the matrices are real. For instance, if a < and < b < c then the first 
two matrices are real. In that case, the cubic is a Helton- Vinnikov curve, and its 
bounded region, when drawn in the affine plane {x = 1}, is the spectrahedron 



1 



ay 



1 + by 







l + cy 



y o 



The symbol 'V" means that the matrix is positive semidefinite. This spectrahed 
is depicted in Figure [T] for the parameter values a = — 1, b = 1 and c — 2. 
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Figure 1. A cubic Helton- Vinnikov curve and its spectrahedron 

This article is organized as follows. In Section 2 we translate into a system 
of polynomial equations in the matrix entries of A,B,C, we determine the num- 



ber of solutions (in Theorem 2.1), and we discuss practical aspects of computing 
these solutions using both symbolic and numeric software. Section 3 is devoted to 
geometric constructions for obtaining the representation ([!]) . Following Dixon [10] , 
these require finding contact curves of degree d — 1 for the given curve of degree d. 

An explicit formula for appears in the article of Helton and Vinnikov 
[161 Eq. 4.2]. That formula requires the numerical evaluation of Abelian integrals 
and theta functions. In Section 4, we explain the Helton- Vinnikov formula, and 
we report on our computational experience with the implementations of [TJ [HI IH] 
in the Maple package algcurves. In Section 5 we focus on the case of quartic 
polynomials and relate our results in |21j to the combinatorics of theta functions. 
Smooth quartics have 36 inequivalent representations ([I]). In the Helton- Vinnikov 
case, twelve of these are real, but only eight are real definite. One of our findings 
is an explicit quartic in Q[x, y, z) that has all of its 12 real representations over Q. 



2. Solving Polynomial Equations 

Our given input is a homogeneous polynomial f(x, y, z) of degree d, usually over Q. 
We assume for simplicity that the corresponding curve in the complex projective 
plane is smooth, we normalize so that f(x, 0, 0) = x d , and we further assume that 
the factors of the binary form f(x,y,0) — Yl i=1 (x + fay) are distinct. Under these 
hypotheses, every equivalence class of representations ([I]) contains a representative 
where A is the identity matrix and B is the diagonal matrix with entries f3\ < 02 < 
• • • < fid . This follows from the linear algebra fact that any two quadratic forms 
with distinct eigenvalues can be diagonalized simultaneously over C. See Section 
IX. 3 in Greub's text book [T3] or the proof of Theorem 4.3 in [2"T] . 
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After fixing the choices A = diag(l, 1, . . . , 1) and B — diag(/?i, 02, •••■> (3d) for 
the first two matrices, we are left with the problem of finding the ( entries of 
the symmetric matrix C — (c«j). By equating the coefficients of all terms x a y^ z 1 
with 7 > 1 on both sides of dlb, we obtain a system of (f^ 1 ) polynomial equations 
in the ( d ^ 1 ) unknowns cy. More precisely, the coefficient of x a y 13 ' z 1 in dlb leads to 
an equation of degree 7 in the cy . We are thus faced with the problem of solving a 
square system of polynomial equations. The expected number of complex solutions 
of that system is, according to Bezout's Theorem, 

x d . 2 d-l . 3 d-2 . 4 d-3 (d-lf-d. (2) 

This estimate overcounts the number of equivalence classes of representations ([!]) 
because we can conjugate the matrix Ax + By + Cz by a diagonal matrix whose 
entries are +1 or —1. This conjugation does not change A or B but it leads to 2 d_1 
distinct matrices C all of which are equivalent. Hence, we can expect the number 
of inequivalent linear matrix representations ([I]) to be bounded above by 

3 d - 2 • 4 d - 3 (d - l) 2 • d. (3) 

We shall refer to this number as the Bezout bound for our problem. 

It is a result in classical algebraic geometry that the number of complex 
solutions to our equations is finite, and the precise number of solutions is in fact 
known as well. The following theorem summarizes both what is known for arbitrary 
smooth curves over C and what can be shown for Helton- Vinnikov curves over R: 

Theorem 2.1. The number of equivalence classes of linear symmetric determinantal 
representations ^ of a generic smooth curve of degree d in the projective plane is 

2 (V)-i.( 2 (V) + i), (4 ) 

unless d > 11 and d is congruent to ±3 modulo 8, when the number drops by one. 
In the case of a Helton- Vinnikov curve, the number of real equivalence classes of 
symmetric linear determinantal representations M is either 
or one less. The number of real definite equivalence classes is precisely 2v 2 /. 

Sketch of Proof. The equivalence classes of representations correspond to in- 
effective even theta characteristics [3] on a smooth curve of genus g = ( d ^ 1 )- The 
number of even theta characteristics is 2 9 ~ 1 (2 9 + 1), and all even theta character- 
istics are ineffective for d < 5 and d = 0, 1, 2, 4, 6, 7 mod 8. In all other cases there 
is precisely one effective even theta characteristic, provided the curve is generic. 
This was shown by Meyer-Brandis in his 1998 diploma thesis [19], and it refines 
results known classically in algebraic geometry [TTJ Chapters 4-5]. The count of 
real and real definite representations will be proved at the end of Section [4] □ 



Computing Linear Matrix Representations 



■5 



The following table lists the numbers in ^ and Q for small values of d: 

degree d 2 3 4 5 6 7 

genus g 1 3 6 10 15 

Bezout bound 1 3 36 2160 777600 1959552000 

True number 1 3 36 2080 524800 536887296 

This table shows that computing all solutions to our equations is a challenge when 
d > 6. Below we shall discuss some computer experiments we conducted for d < 5. 

As before, we fix A to be the dxd identity matrix, denoted Id^, and we fix 
B to be the diagonal matrix with entries fix < • • • < fid. We also fix the diagonal 
entries of C since these are determined by solving the d linear equations that arise 
by comparing the coefficient of any of the d monomials x l y d ~ i ~ l z in Q. They are 
expressed in terms of / and the ft by the following explicit formula: 



|£(-a,i,o) 
% (-ft, 1,0) 



ft- a ') „ ' far» = l,2,...,d. (5) 



We are thus left with a system of („) equations in the (f\ off-diagonal un- 
knowns Cij. In order to remove the extraneous factor of 2 d ~ 1 in the Bezout bound 
([2]) coming from sign changes on the rows and columns of C, we can perform a 
multiplicative change of coordinates as follows: Xij — cfj for j — 1,2, ... ,d and 
Cij for 2 < i < j < d. This translates our system of polynomial equa- 
tions in the c^- into a system of Laurent polynomial equations in the i, 3 , and each 
solution to the latter encodes an equivalence class of 2 d_1 solutions to the former. 

Example 2.2. Let d = 4. We shall illustrate the two distinct formulations of the 
system of equations to be solved. We fix a quartic Helton- Vinnikov polynomial 



f{x,y,z) = dot 



x + 0\y + 7112 7i2z 713Z 7142 

7122 x + /3 2 y + 7222 7232 7242 

7132 7232 x + /3 3 y + 7332 7342 

7142 7242 7342 x + 4 y + 7442 



where and jjk are rational numbers. From the quartic f(x,y,z) alone we can 
recover the ft and the diagonal entries jjj as described above. Our aim is now to 
compute all points (012,013,014,023,024,034) € C 6 that satisfy the identity 



dct 



x + Piy + fnz C122 C132 C142 

C122 x + P2JI + 7222 C232 C242 

C132 C232 x + fay + J33Z c 3i z 

C142 C242 C342 x + j3iy + 7442 



= f(x,y,z). 



The coefficient of z gives one equation of degree 4 in the six unknowns Cjj , the co- 
efficients of xz 3 and yz 3 give two cubic equations, and the coefficient of x 2 z 2 , xyz 2 
and y 2 z 2 give three quadratic equations in the cy. The number of solutions in C 6 
to this system of equations is equal to 2 3 3 2 4 = 288. These solutions can be found 
using symbolic software, such as Singular [6]. However, the above formulation 
has the disadvantage that each equivalence class of solutions appears eight times. 
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We note that, for generic choices of j3i, 7jfc, all solutions lie in the torus (C*) 6 
where C* = C\{0}, and we shall now assume that this is the case. Then the 8-fold 
redundancy can be removed by working with the following invariant coordinates: 

x 12 = c\ 2 , X13 = C 2 3 , X14 = C 2 4 , 

^23 = C12C13C23 , X 2i = Ci 2 C 14 C 2 4 , X34 = C13C14C34- 

We rewrite our six equations in these coordinates by performing the substitution: 



C12 



C23 



x 



1/2 

12 j c 13 

X23 

1/2 1/2 ' 
X 1 X 1 



X 



1/2 
13 ■ 



C 2 4 



C14 = X 

%24 
1/2 1/2 
X 12 X 14 



1/2 
14 ' 



C34 



£34 



1/2 1/2' 
x 13 X 14 



b 12 ^13 

This gives six Laurent polynomial equations in six unknowns x 12l £13, X14, X23, X24, 
X34. They have precisely 36 solutions in (C*) 6 , one for each equivalence class. □ 

While the solution of the above equations using symbolic Grobner-based soft- 
ware is easy for d = 4, we found that this is no longer the case for d > 5. For d = 5, 
it was necessary to employ tools from numerical algebraic geometry, and we found 
that Bertini [3] works well for our purpose. The computation reported below is 
due to Charles Chen, an undergraduate student at UC Berkeley. This was part of 
Chen's term project in convex algebraic geometry during Fall 2010. 

For a concrete example, let us consider the following polynomial which defines 
a smooth Helton- Vinnikov curve of degree d = 5: 

f(x, y, z) = x 5 + 3x 4 y - 2x 4 z - 5x 3 y 2 - I2x i z 2 - 15x 2 y 3 + Wx 2 y 2 z - 28x 2 yz 2 + 14x 2 z 3 + 
Axy 4 - 6xy 2 z 2 - \2xyz 3 + 26zz 4 + 12j/ 5 - 8y i z - 32y :i z 2 + 16y 2 z 

The symmetric linear determinantal representation we seek has the form 




~ 3 +48i/z 4 -24z 5 



x + y 







- 2;/ 



2y 








3y- 



2z 



+ C-z, 



where C — (cy) is an unknown symmetric 5x5-matrix with zeros on the diagonal. 
This leads to a system of 10 polynomial equations in the 10 unknowns c^-, namely, 
4 quadrics, 3 cubics, 2 quartics and one quintic. The number of complex solutions 
equals 16 • 2080 = 33280, which is less than the Bezout bound of 2 4 • 3 3 • 4 2 • 5 = 
16-2160 = 34560. One of the 33280 solutions is the following integer matrix, which 
we had used to construct f(x, y, z) in the first place: 



C 



Of course, the other 15 matrices in the same equivalence class have the same 
friendly integer entries. The other 16 • 2079 = 33264 complex solutions were found 






2 


1 








2 











1 


1 








2 


1 








2 





-1 





1 


1 


-1 
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numerically using the software Bertini [4]. Of these, 16 • 63 are real. Chen's code, 
based on Bertini, outputs one representative per class. One of the real solutions is 





1.8771213868 
0.1333876113 
0.3369345269 
0.2151885297 



1.8771213868 


1.3262201851 
0.1725327846 
1.0570303927 



0.1333876113 
1.3262201851 


-2.0093203944 
-0.8767796987 



0.3369345269 
0.1725327846 
-2.0093203944 


-0.7659896773 



0.2151885297 ' 
1.0570303927 
-0.8767796987 
-0.7659896773 




3. Constructing Contact Curves 



A result in classical algebraic geometry states that the equivalence classes of sym- 
metric linear determinantal representations of a plane curve of degree d are in 
one-to-one correspondence with certain systems of contact curves of degree d— 1. 
Following IT] Prop. 4.1.6], we now state this in precise terms. Suppose that our 
given polynomial is / = det(-M) where M = is a symmetric dxd- matrix of 
linear forms in x, y, z. We can then form the dxd adjoint matrix, adj(M), whose 
entry my is the (z, j)th (d — l)-minor of M multiplied by (— iy +J . For any vector 
of parameters u = (ui, U2, ■ ■ . , Ud) T , we consider the degree d — 1 polynomial 



g u (x,y,z) 





Ul 


T r 








u T ad]{M)u = 








Ud_ 





m i2 



m i2 
m 22 



m d i m d2 



mid 
m 2 d 

mdd 





Ul 








Ud_ 



(6) 



The curve V{g u ) has degree d—1, and it is a contact curve, which means that all in- 
tersection points of V(f) and V(g u ) have even multiplicity, generically multiplicity 
2. To see this, we use [HJ Lemma 4.1.7], which states that, for any u,v E C d , 

g u (x,y,z) ■ g v {x,y,z) - (u T adj(M» 2 G (/) in C[x,y,z]. 

In particular, for u = e^, v = ej, this shows that both V(ma) and V{rrijj) are 
contact curves, and V(my) meets V(/) in their d{d— 1) contact points. 

We say that two contact curves V(gi) and V(<72) of degree r lie in the same 
system if there exists another curve V(h) of degree r that meets V(/) precisely in 
the r-d points V(/, gi) U V(/, g 2 )- A system of contact curves is called syzygetic if 
it contains a polynomial of the form £ 2 g, where £ is linear and g is a contact curve 
of degree r — 2, and azygetic otherwise. A contact curve of V(/) is called syzygetic, 
resp. azygetic, if it lies in a system that is syzygetic, resp. azygetic. 

Dixon [TO] proved that the contact curves V(g u ) are azygetic and all azygetic 
contact curves of degree d—1 appear as g u for some determinantal representation 
/ = det(M). In particular, he gives a method of constructing a determinantal 
representation M for / starting from one azygetic contact curve of degree d—1. 

The input to Dixon's algorithm is an azygetic contact curve g of degree d—1 
of the given curve / of degree d. Given the two polynomials / and g, the algorithm 
constructs the matrix M — adj(M) in It proceeds as follows. Since V{g) meets 
V(/) in d(d— l)/2 points, the vector space of polynomials of degree d—1 vanishing 
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Figure 2. A quartic Helton- Vinnikov curve and four contact cubics 



at these points (without multiplicity) has dimension d. Let Wn = g and extend 
mn to a basis {mn,mi2, . . . , mid} of this vector space. For i,j G {2, 3, . . . ,d}, 
the polynomial mumij vanishes to order two on V(/,mn), so it lies in the ideal 
(mix,/). Using the Extended Buchberger Algorithm, one finds a degree d — 1 
polynomial m^- such that mumy — mumij G (/). The cZxcZ-matrix M = (m^) 
has rank 1 modulo (/), therefore its 2x2-minors are multiples of /. This implies 
that the adjoint matrix adj(M) = det(M) • M" 1 has the form \f d ~ 2 ■ M where 
A G C\{0} and M is a symmetric matrix of linear forms with det(M) = /. One 
could run through this construction starting from a syzygetic contact curve, but 
the resulting matrix M would have determinant zero. 

The main challenge with Dixon's algorithm is to construct its input polyno- 
mial g. Suitable contact curves are not easy to find. A symbolic implementation of 
the algorithm may involve large field extensions, and we found it equally difficult 
to implement numerically. For further discussions see [IH1 §2.2] and [HJ §2]. 

Remark 3.1. Starting from a real azygetic contact curve g, one can use Dixon's 
method to produce a real determinantal representation of /. A determinantal 
representation M is equivalent to its conjugate M if and only if the system of 
contact curves {w T adj(Af)u : u G C 4 } C C[x,y, z\% is real, i.e. invariant under 
conjugation. The representation M is equivalent to a real matrix if and only if this 
system contains a real contact curve. By [14, Prop 2.2], if the curve V(/) has real 
points, then these two notions of reality agree. However, this approach does not 
easily reveal whether an equivalence class contains a real definite representative. 

Example 3.2. The following Helton- Vinnikov quartic was studied in [21] Ex. 4.1]: 



r,4, 4 , 4 o22 o22, 22 

= 2x +y + z — ox y — ox z + y z . 
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It is shown on the left in Figure [2] It has a symmetric determinantal representation 



f{x,V,z) 



dot 



ux + y 





az 


bz 





ux — y 


cz 


dz 


az 


cz 


x + y 





bz 


dz 





x — 



(7) 



where a = -0.5746..., b = 1.0349..., c = 0.6997..., d = 0.4800... and u = y/2 are 

the coordinates of a real zero of the following maximal ideal in Q[a, b, c, d, u]: 

(u 2 - 2,256d 8 - 384d 6 u+256d 6 -384d 4 u+672d 4 -336d 2 u+448d 2 -84u +121, 



23c- 
23&- 



7584d 7 
5760d 7 i 



T0688d 7 -5872d 5 u-8384d 5 +1806d 3 w+2452d 3 -181dit-307d, 
-8192d 7 -4688d 5 u-6512d 5 -H452d 3 u+2200d 3 -212dM-232d, 



Hes(/) = det 



23a - 1440d 7 M -2048d 7 +1632d 5 u+2272d 5 -570d 3 u-872d 3 +99du+81d). 

The principal 3x3-minors of the 4x4-matrix in Q are Helton- Vinnikov polynomi- 
als of degree 3. They are the four contact cubics shown on the right in Figure[2j □ 

In summary, Dixon's method furnishes an explicit bijection between equiva- 
lence classes of symmetric determinantal representations (JlJ of a fixed curve V(/) 
of degree d and azygetic systems of contact curves of V(/J of degree d — 1. 

For d = 3, there is another geometric approach to finding representations (JlJ. 
We learned this from Didier Henrion who attributes it to Frederic Han. Suppose 
we are given a general homogeneous cubic /(a;, y, z). We first compute the Hessian 

d 2 f/dx 2 d 2 f/dxdy d 2 f/dxdz 
d 2 f/dxdy d 2 f/dy 2 d 2 f/dydz 
d 2 f/dxdz d 2 f/dydz d 2 f/dz 2 

This is also a cubic polynomial, and hence so is the linear combination i-/+Hes(/), 
where t is a parameter. We now take the Hessian of that new cubic, with the aim 
of recovering /. It turns out that we can do this by solving a cubic equation in t. 

Proposition 3.3. There exist precisely three points (s,t) 6 C 2 such that 

s-f = Hes(i-/ + Hes(/)). (8) 
The resulting three symmetric determinantal representations of f are inequivalent. 

Proof. The statement is invariant under linear changes of coordinates in P 2 , so, 
by [2 Lemma 1] , we may assume that the given cubic is in Hesse normal form: 

f(x,y,z) = x 3 + y 3 + z 3 — mxyz. 

In that case, the result follows from the discussion in fL7\ page 139]. Alternatively, 
we can solve the equations obtained by comparing coefficients in (|8|). This leads to 



(12m 4 + 2592m)t - 16m 6 + 8640m 3 + 93312 = 
(12m 4 + 2592m)t 2 + (48m 6 - 25920m 3 - 279936)4 



and 

48m 8 + 20736m 5 



- 2239488m- 3 



This has three solutions (s, t) g C 3 . The resulting representations ([l]) are inequiva- 
lent because the Hessian normal form of a PGL(3, C)-orbit of cubics is unique. □ 
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4. Evaluating Theta Functions 

The proof of the Helton- Vinnikov Theorem relies on a formula, stated in [Tr?l 
Eq. 4.2], that gives a positive definite determinantal representation of a Helton- 
Vinnikov curve in terms of theta functions and the period matrix of the curve. 
Our aim in this section is to explain this formula and to report on computational 
experiments with it. Numerical algorithms for computing theta functions, period 
matrices and Abelian integrals have become available in recent years through work 
of Bobenko, Deconinck, Heil, van Hoeij, Patterson, Schmies, and others [3 [HI E] - 
There exists an implementation in Maple, and we used that for our computations. 
Our Maple worksheet that evaluates the Helton- Vinnikov formula can be found at 

www . math . uni-konstanz . de/ ^plaumann/ theta . html (9) 

Before stating the Helton- Vinnikov formula, we review the basics on theta 
functions. Our emphasis will be on clearly defining the ingredients of the formula 
rather than explaining the underlying theory. For general background see |20j . Fix 
jeN and let H 9 be the Siegel upper half-space, which consists of all complex, sym- 
metric g x (^-matrices whose imaginary part is positive definite. The Riemann theta 
function is the holomorphic function on C 9 x % 9 defined by the exponential series 

9{\i, Q) = exp(7ri(m T f2m + 2m T u)) , 

where i = y— T, u = (ui,...,u g ) € C 9 and fl € T-L 9 . We will only need to 
consider 9(u, Q) as a function in u, for a fixed matrix f2, so we may drop Q from 
the notation. In other words, we define 9: C 9 — >• C by 0(u) = 9(u, SI). The theta 
function is quasi-periodic with respect to the lattice 7L 9 + ilZ 9 C C 9 , which means 

6»(u + m + £!n) = exp(7ri(-2n T u - n T fln)) ■ 9(u), for all m, n E Z 9 . 

A theta characteristic is a vector e = a+fib e C 9 with a, b g {0, |} 9 . The function 

6>[e](u) = exp(vri(a T r2a + 2a T (u + b))) -6l(u + f7a + b) 

is the theta function with characteristic e. There are 2 29 different theta character- 
istics, indexed by ordered pairs (2a, 2b) of binary vectors in {0, l} 9 . We also use 
T2al 

the notation 9 ^ (u) for the function 9[e](u). For e = we simply recover 9(u). 

Let / <E C[x,y,z]d be a homogeneous polynomial of degree d. Assume that 
the projective curve X = Vc(/) is smooth and thus a compact Riemann surface of 
genus g = \{d— l)(e? — 2). Let (ui^, . . . ,u>' g ) be a basis of the ^-dimensional complex 
vector space of holomorphic 1-forms on X , and let a±, . . . , a g , Pi, . . . , /3 g be closed 
1-cycles on X that form a symplectic basis of Hi(X,Z) = Z 9 x Z 9 . This means 
that the intersection numbers of these cycles on X satisfy (a.j ■ ctk) = {Pj • Pk) = 0, 
and (ctj ■ Pk) = 8jk for all j, k = 1, . . . , g. The period matrix of the curve X with 
respect to these bases is the complex g x 2g-matrix {VL\ \ f^) whose entries are 

(fli)jk = / u'j and (fi 2 )jfc = / , for j, k = 1, . . . , g. 



Computing Linear Matrix Representations 



11 



The g x ^-matrices fii and 0,% are invertible. Performing the coordinate change 
(wi,...,w s ) = (u[,...,L)' g ) ■ (fl^ 1 ) 7 

leads to a basis in which the period matrix is of the form ( ld g | fi] -1 ^ ) • The basis 
lu = (lui, . . . , ujg) is called a normalized basis of differentials and depends uniquely 
on the symplectic homology basis. The g x g-matrix i7 = r2^ 1 J7 2 is symmetric and 
lies in the Siegel upper half-space % 9 . It is called the Riemann period matrix of the 
polynomial f(x, y, z) with respect to the homology basis (a\, . . . , a g ,/3i, . . . , f3 g ). 

With the given polynomial / we have now associated a system {#[e]( • , O)} of 
2 2g theta functions with characteristics. A theta characteristic e = a + ilh is called 
even (resp. odd) if the scalar product (2a) T (2b) of its binary vector labels is an even 
(resp. odd) integer. This is equivalent to 9[e] being an even (resp. odd) function in 
u. In symbols, we have #[e](— u) = ( — l) 4a b 6>[e](u). Changing symplectic bases 
of Hi(X,Z) corresponds to the right-action of the symplectic group Sp 2ff (Z) on 
the period matrix (Qi | 2 ). This action will permute the theta characteristics. In 
particular, there is no distinguished even theta characteristic 0. 

Finally, we define the Abel-Jacobi map by <fi(P) = (J Po uii, ■ ■ ■ ,J Pq u g ) T for 
Pel, where Pq G X is any fixed base point. This is a holomorphic map, but it 
is well-defined only up to the period lattice A = Z 9 + SXL 9 C C 9 . In other words, 
the Abel-Jacobi map is a holomorphic map tj>: X — > Jac(A) = C 9 /A. 

We are now ready to state the formula for ([lj in terms of theta functions. 

Theorem 4.1 (Helton- Vinnikov [16 ). Let f € R[x, y, z] d with /(l, 0, 0) = 1 and let 

X denote Vfc(/) C P 2 . We make the following two assumptions: 

1. The curve X is a non-rational Helton-Vinnikov curve with the point (1:0:0) 



e 



inside its innermost oval. The latter means that, for all v € M 3 \{0}, th 
univariate polynomial f(v + t ■ (1, 0, 0)) G R[t] has only real zeros. 
2. The d real intersection points of X with the line {z = 0} are distinct non- 
singular points Qi, . . . , Qd, with coordinates Qi = (— f3j : 1 : 0) where j3j =/= 0. 

Then f(x, y, z) = dct{\AdX + By + Cz) where B — diag(/3i, . . . , (3d) and C is real 
symmetric with diagonal entries Cjj as in The off-diagonal entries of C are 



p k -fj eis^jQ,)-^)) o il Vg[gK0) \ u-Ve\e](Q) 

Cjk ~ e[W)'m{<KQ*)-<KQj))\ ~ d ^/v) J \ - d ( z /y) 

Here e is an arbitrary odd theta characteristic and 5 is a suitable even theta charac- 
teristic with 9[8\{Q) 7^ 0. The theta functions are taken with respect to a normalized 
basis of differentials uj — (oji, . . . , u> g ), and <f>: X — ¥ Jac(A) is the Abel-Jacobi map. 



The remarkable expression for the constants Cjk in ( 10 1 does not depend on 
the choice of the odd characteristic e. If the curve X is smooth, then all equiv- 
alence classes of symmetric detcrminantal representations are obtained when 5 
runs through all non-vanishing even theta characteristics. The proof of Theorem 
4.1 given in [16] is only an outline. It relies heavily on earlier results on Riemann 
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surfaces due to Ball and Vinnikov in [2, 24]. As we found these not easy to read, we 
were particularly pleased to be able to verify Theorem |4 . 1 1 with our experiments. 



The Helton- Vinnikov formula ( 10 ) remains valid when X is a singular curve. 
In that case the period matrix, the differentials, and the Abel-Jacobi map are 
meant to be defined on the desingularization of A, a compact Riemann surface of 
genus g with g < h(d — l)(d — 2). The formula holds as stated, but one no longer 
obtains all equivalence classes of symmetric determinantal representations. 

The Riemann period matrix, the theta functions, their directional derivatives, 
and the Abel-Jacobi-map can all be evaluated numerically in recent versions of 
Maple. When computing the expressions under the square roots, note that both 
the numerator and denominator are 1-forms on X. Every holomorphic 1-form 
on the curve X can be written as r ■ du, where u = z/x, v = y/x, and r is a 
rational function in u and v. The algcurves package in Maple will compute uj in 
this form, so we obtain u)j — rj(u,v) ■ du. To evaluate the 1-form d(z/y), we set 

^du+^( 

ou ov 



h(u, v) = f(l,v, u) and use the identity dh(u, v) — ^du + j^dv — 0. This implies 



I v ( 1 u— 

d(z/y) = d(u/v) = -du-^dv = [-+to] d (n) 

v ov 



so that d(z/y)(Qj) = —fijdu. Under the square root signs in (10 1, the factor du 



appears in the numerator u> and also in the denominator (111, and we cancel 



it. Hence the expressions under the square roots are rational functions, namely 



r(u,v) ■ V0[e](O) divided by the expression in parentheses on the right in (11), 
where r(u, v) is the vector of rational functions Tj(u, v). 

While the evaluation of theta functions is numerically stable, we found the 
computation of the period matrix and the Abel-Jacobi map to be more frag- 
ile. Computing the d vectors 4>{Qj) is a l so by far the most time-consuming step. 
Nonetheless, Maple succeeded in correctly evaluating the Helton- Vinnikov formula 
for a wide range of curves with d < 4, and for some of degree d — 5. However, the 



off-diagonal entries in ( 10 1 we found in our computations were sometimes wrong by 
a constant factor (independent of j, k), for reasons we do not currently understand. 
For a concrete example take the quartic in Example |3.2| Using the formula 



( 10 1 we obtained all eight definite determinantal representations detflddir. + By + 
Cz). Our Maple code runs for a few minutes and finds all solutions accurately with 
a precision of 20 digits. We verified this using the prime ideal in Example |3.2| 

The case of smooth quintics (genus 6) is already a challenge. With the help of 
Bernard Deconinck, we were able to compute a determinantal representation (jlOh 
with an error of less than 10 -3 for the quintic polynomial at the end of Section [2l 



However, the representation we obtained was not real (see Remark 4.4 1. 

We conclude this section with the proof of the second part of Thm. 2.1 and 
discuss what are the suitable choices for the even theta characteristic S in Thm. l4~T1 
that will lead to real and real definite equivalence classes of representations. Note 
that a representation obtained from the theorem is real definite if and only if 
the matrix C is real. The real non-definite equivalence classes of representations 
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correspond to the case when C is a non-real matrix for which there exists a matrix 
U e GL d {C) such that UU T , UBU T , and UCU T have all real entries. Whether 
such U exists for given complex symmetric matrices B and C is not at all obvious. 
An explicit example is given in Ex. |5.2| 



For the proof of Thm. 2.1 we repeat the relevant part of the statement: 

Theorem 4.2. Let f € R[x, y, z]d be homogeneous and assume that the projective 
curve X — Vc(/) is a smooth Helton-Vinnikov curve. The number of real equiva- 
lence classes of symmetric determinantal representations is generically either 

2 9- 1 (2 fc - 1 + 1) 

or one less, where k = |~|] is the number of connected components of the set of 
real points X(R). Of these real equivalence classes, exactly 2 9 are definite. 

Proof. The result follows from work of Vinnikov |24j on self-adjoint determinantal 
representations, which we apply here to our situation. By [2H Prop. 2.2], a sym- 
plectic basis of H\(X, Z) can be chosen in such a way that the Riemann period 
matrix f2 satisfies Q + ft = H, where H is a g x g block diagonal matrix of rank 

in the top left corner and all other entries zero. 



g—k+1 with r/2 blocks 



1 



4.1 



The linear symmetric determinantal representation obtained by Thm 
from an even theta characteristic 5 is equivalent to a real one if and only if S is 
real, i.e. invariant under the action of complex conjugation on the g-dimensional 
torus Jac(A). Since X(M.) ^ 0, any conjugation-invariant divisor class on X con- 
tains a real divisor (see [121 Prop. 2.2]). From such a divisor, one can construct 
a symmetric determinantal representation (see [5] or [S3] ). When the symplectic 
basis of Hi(X, Z) is chosen as above, the action of complex conjugation on Ja,c(X) 
is given byC = u + f2vi-^( : =u + Q(Hu - v) (see [231 Prop. 2.2]). For the 
even theta characteristic 6 = a + fib, a, b g {0, |} 9 , the condition 5 — S in 
Jac(X) = C 9 /(Zf + SO?) thus becomes 

Ha = 2b mod Z 5 . 

This happens if and only if a% = ■ ■ ■ = a r = 0. Counting the possible choices 
of a r+ i, . . . , a g and b, we conclude that there are exactly 2 9-1 (2 fe ~ 1 + 1) even 
real theta characteristics. From the first part of Thm. |2.1[ we know that when 
d = ±3 mod 8, exactly one even theta characteristic vanishes, i.e. 0[<$](O) = 0. 
All other even theta characteristics are non-vanishing and therefore correspond 
to determinantal representations. Furthermore, by |24[ Thm. 6.1], an even theta 
characteristic 6 = a + fib will correspond to a real definite equivalence class if and 
only if a = 0, and all such 5 are always non- vanishing [2H Cor. 4.3]. Thus there 
are exactly 2 9 definite representations, since b can be any element of {0, . □ 

Example 4.3. When g = 3 and a homology basis has been picked as above, the 
even real theta characteristics are given by the 12 binary labels 



000 




000 




000 




000 




000 




000 




000 




000 




001 




001 




001 




001 


000 




001 




010 




011 




100 




101 




110 




111 




000 




010 




100 




110 



(12) 
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The first eight labels correspond to the definite representations and the last four 
correspond to the non-definite real equivalence classes of representations. 

Remark 4.4. The characterization of real and real definite even theta character- 
istics provided by the proof of Thm. |4.2| depends on the choice of a particular 
symplectic homology basis. Unfortunately, the current Maple code for computing 
the period matrix does not give the user any control over the homology basis. This 
makes it hard to find real representations using Thm. |4.1| in any systematic way. 



5. Quartic Curves Revisited 

In this section we focus on the case of smooth quartic curves, studied in detail 
in |21j . so we now fix d = 4 and g — 3. Quartic curves are special because they 
have contact lines, i.e. bitangents, and we can explicitly write down higher degree 
contact curves as products of bitangents. This was exploited in [5TJ §2], where we 
used ayzgetic triples of bitangents as our input to Dixon's algorithm (Section 3). 

Plane quartics are canonical embeddings of genus 3 curves |llj . and there is 
a close connection between contact curves and theta functions. The 28 bitangents 
of the curve are in bijection with the 28 odd theta characteristics e = a + Qh, and 



this will be made explicit in ( 15 ) below. The 36 azygetic systems of contact cubics 
correspond to the 36 even theta characteristics. As seen in Example |4.3| of the 
resulting 36 determinantal representations, 12 are real, but only 8 are definite. We 
can also derive the number 12 from the combinatorics of the bitangents as in |21j . 

Proposition 5.1. A smooth Helton- Vinnikov quartic V(/) has exactly 12 inequiva- 
lent representations f = det(Az; + By + Cz) with A, B, C symmetric and real. 



Proof. This is a special case of Theorem |2.1[ however, we here give an alternative 
proof using the setup of |21j . Let M be a symmetric linear determinantal represen- 
tation of / and M. the system of contact cubics {it T adj(Af )u} C C[x,y,z]3. The 
representation M is equivalent to its conjugate M if and only if the system M. is 
real, i.e. invariant under conjugation. The representation M is equivalent to a real 
matrix if and only if M. contains a real cubic. The matrix M induces a labeling 
of the 28 = (2) bitangents, bij, with 1 < i < j < 8. The system M. is real if and 
only if conjugation acts on the bitangents via this labeling, that is, there exists 
7r G S"g such that bij = b n u\ n u\. Since / is a Helton- Vinnikov polynomial, this 
permutation will be the product of four disjoint transpositions (see |2TJ Table 1]). 

Suppose A4 is real, with permutation tt € S%. The other 35 representations ([I]) 
correspond to the m /2 partitions of {1, . . . , 8} into two sets of size 4. If I\I C is such 
a partition then the corresponding system of contact cubics contains 56 products 
of three bitangents, namely bijbikbu and bi m bj rn bk£ where i,j,k,l,m are distinct 
and {i,j,k,l} = I or I c . This system is real if and only if tt fixes the partition 
I\I C . There are exactly 11 such partitions: if it = (12)(34)(56)(78), they are 

1234|5678, 1256|3478, 1278|3456, 1357|2468, 1358|2467, 1368|2457 , , 
1367|2458, 1457|2368, 1458|2367, 1467|2358, and 1468|2357. [ ' 
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Together with the system M., there are 12 real systems of azygetic contact cubics. 

Next, we will show that each of these systems actually contains a real cubic. 
To do this, we use contact conies, as the product of a bitangent with a contact conic 
is a contact cubic. By [21] Lemma 6.7], there exists a real bitangent b <= M[x, y, z]i 
and a real system of contact conies Q C C[x, y, z]2 such that their product {b ■ q : 
q € Q} lies in the system M C C[x,y, z] 3 . Furthermore, by [HJ Prop. 6.6], since 
Vk(/) is nonempty, every real system of contact conies Q to / contains a real conic 
q. The desired real contact cubic is the product b ■ q. □ 



The technique in the last paragraph of the above proof led us to the following 
result: There exists a smooth Helton- Vinnikov quartic / € Q[x,y, z]^ that has 12 
inequivalent matrix representations ([!]) over the field Q of rational numbers. 

Example 5.2. The special rational Helton- Vinnikov quartic we found is 

f(x,y,z) = 93081a; 4 + 53516ir 3 y - 73684zV + -31504a;y 3 + 9216y 4 

- 369150a; 2 z 2 - 15970Cte?/z 2 + 57600y 2 z 2 + 90000z 4 . 

This polynomial satisfies f(x,y,z) = det(-M) where 



M = 



50a; 
-25a; 
-26a; - 34y - 25z 
9a; + 6y + 15z 



-25a; 
25a; 

27x + 18y - 20z 
—9x — 6y 



-26x - 34t/ - 25z 
27a; + 18y - 20z 
108a; + 72y 
-18a; - 12y 



9x + 6y + 15z 
—9x — 6y 
-18a; - 12y 
6x + 4:y 



This representation is definite because the matrix M is positive definite at the 
point (1:0: 0). Hence V(/) is a Helton- Vinnikov curve with this point in its inner 
convex oval. Rational representatives for the other seven definite classes are found 
at our website ([9]), along with representatives for the four non-definite real classes. 
One of them is the matrix 



M 



1 KiS 



25x 


-32x + 12y 
-60z 





25a; 
lOz 
24a; + 16y 



-32a; + I2y 
lOz 
6a; + 4y 




-60z 
24a; + 16y 


6x + Ay 



(14) 



We have det(M!4g 8 ) = 4 • f{x, y, z), and this matrix is neither positive definite nor 
negative definite for any real values of z. Any equivalent representation of a 
multiple of / in the form det^diX+By+Cz) considered in Sections [2] and [1] cannot 
have all entries of C real. One such representation, for a suitable U £ GL^C), is 



U'MuesU = 



7i y 







- T || ?j v / 26980 ^ 
t^II V16330 z 



-^1^/570 2 



T ||gV / 26980 i 2 

o 



^^16330 2- 
4v / 345 i 2 



2:i 







23 y 



The correspondence between bitangents and odd theta characteristics can be 
understood abstractly via the isomorphism between the Jacobian of X and the 
divisor class group Gl°(X), and can be turned into an explicit formula for the 
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bitangents. Let u — z/x, v — y/x and write h(u,v) = f(l,v,u). Then a basis for 
the 3-dimensional complex vector space of holomorphic 1-forms on X is given by 

/ du vdu udu 



\dh/dv' 8h/dv> dh/dv 

Let (Oil ) be the period matrix of / with respect to (lji, u>2, W3) and any sym- 
plectic basis of H 1 (X,Z). Then given an odd theta characteristic e = a + fib, the 
corresponding bitangent is defined by the linear form 



b e {x,y,z) 



(x,y,z) 



(15) 



(V%](0)) J 

where V#[e] is the gradient of 9[e] in the three complex variables z\,Z2,z$. For the 
proof, see Dolgachev [HI Section 5.5.4]. This holds independently of the symplectic 
basis of Hi(X, Z), but a change of that basis will permute the bitangents. 



The formula (15 1 can be evaluated using the Maple code described in Sec- 
tion 4. This allows us to compute the 8 x 8-bitangent matrix (fry) of [2TJ Eq. 3.4] 
directly from the Riemann period matrix £1 of the curve X, using a technique due 
to Riemann described in [151 §2]. In this manner, one computes the symmetric 
determinantal representations ([lj of the curve X directly from the period matrix 
n. This computation seems to be a key ingredient in constructing explicit three- 
phase solutions of the Kadomtsev-Petviashvili equation [T2] , and we hope that the 
combinatorial tools developed here and in [21] will be useful for integrable systems. 

One of the earliest papers on algorithms for theta functions in genus three was 
written by Arthur Cayley in 1897. In [5] he gives a concrete bijection between the 
bitangents fry of a plane quartic and the odd theta characteristics, and also between 
the classes I\I C of determinantal representations and the even theta characteristics. 
We here reproduce a relabeled version of the table in Cayley's article: 



2b\2a 


000 


100 


010 


110 


001 


101 


011 


111 


000 





1238 


1267 


1245 


1468 


1578 


1356 


1347 


100 


1234 


48 


1235 


35 


1457 


16 


1378 


27 


010 


1256 


1247 


57 


46 


1367 


1345 


23 


18 


110 


1278 


37 


68 


1236 


1358 


25 


14 


1567 


001 


1357 


1346 


1478 


1568 


12 


38 


67 


45 


101 


1368 


26 


1456 


17 


34 


1248 


58 


1235 


011 


1458 


1678 


13 


28 


56 


47 


1257 


1246 


111 


1467 


15 


24 


1348 


78 


1237 


1268 


36 



Here a partition I\I C of {1, . . . , 8} is represented by the 4-tuple / which contains the 
index 1. For instance, the 4-tuple 1238 corresponds to the even theta characteristic 

Each partition I\I C represents a Cremona transformation leading to a new 



100 
000 



representation ([!]) as described in [3TJ §3] . The twelve 4-tuples marked in bold face 
are the real equivalence classes, and this gives a bijection between the lists in ( 12 ) 



and in (13). Likewise, the pairs ij in Cayley's table represent bitangents &y and 



the corresponding odd theta characteristics. For instance, the odd characteristic 
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represents the bitangent 615. In this manner, we can parametrize the 28 
bitangents of all plane quartics explicitly with odd theta functions. 

Experts in moduli of curves will be quick to point out that this parametriza- 
tion should extend from smooth curves to all stable curves. This is indeed the case. 
For instance, four distinct lines form a stable Helton- Vinnikov quartic such as 

f{x,y,z) = xyz(x + y + z). 

The bitangent matrix (bij) of this reducible curve has 7 distinct non-zero entries: 



" 


z 


y 


y + z 


X 


X + z 


x + y 


x+y+z 


z 





y + z 


y 


X + z 


X 


x+y+z 


x + y 


y 


y + z 





z 


x + y 


x+y+z 


x 


X + z 


y + z 


y 


z 





x+y+z 


x + y 


x + z 


X 


X 


X + z 


x + y 


x+y+z 





z 


y 


y + z 


x + 2 


X 


x+y+z 


x + y 


z 





y + z 


y 


x + y 


x+y+z 


X 


X + z 


y 


y + z 





z 


jc+y+z 


x + y 


X + z 


X 


y + z 


y 


z 






All principal 4x4-minors of this 8x8-matrix are multiples of f(x,y,z), most of 
them non-zero. They are all in the same equivalence class, which is real but not 
definite. The entries in the bitangent matrix indicate a partition of the 28 odd 
theta characteristics into seven groups of four. For instance, the antidiagonal entry 
x + y + z corresponds to the four entries 18, 27, 36 and 45 in Cayley's table, and 

hence to the four odd theta characteristics 

If we consider a family of smooth quartics that degenerates to the reducible 
quartic f(x,y,z), then its bitangent matrix will degenerate to the above 8x8- 
matrix, and hence the 28 distinct bitangents of the smooth curve bunch up in 





100 




111 




001 


1 


111 




111 


and 


111 
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seven clusters of four. This degeneration is visualized in Figure [3] Among the 
seven limit bitangents arc the three lines spanned by pairs of intersection points. 

Algebraically, such a degenerating family can be described as a curve over a 
field with a valuation, such as the field of real Puiseux series K{{e}}. The notions of 
spectrahedra and Helton- Vinnikov curves makes perfect sense over the real closed 
field E{{e}}. This has been investigated from the perspective of tropical geometry 
by David Speyer, who proved in |22] that tropicalized Helton- Vinnikov curves 
are precisely honeycomb curves. We believe that the tropicalization in 22J offers 
yet another approach to constructing linear determinantal representations ([I]), in 
addition to the three methods presented here, and we hope to return to this topic. 

Acknowledgments. We wish to thank Charles Chen, Bernard Deconinck, Didicr 
Hcnrion, Chris Swierczewski and Victor Vinnikov for discussions and computa- 
tional contributions that were very helpful to us in the preparation of this article. 
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